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Abstract 



In previous works, we studied the isotropisation of Bianchi 
class A models with a minimally coupled scalar field. In this 
paper, we extend these results to the case of a non minimally 
coupled one. We first make the calculations in the Einstein 
^ I frame where the scalar field is minimally coupled to the cur- 

vature but non minimally coupled to the perfect fiuid. Then, 
I we use a conformal transformation to generalise our results to 

Q^ ' a scalar field non minimally coupled to the curvature. Uni- 

I verse isotropisation for the Brans-Dicke and low energy string 

*0 ' theories are studied. 

o 
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1 Introduction 

X 

^ ' There are numerous reasons to consider the presence of some scalar 

fields in our Universe. Historically, the most famous scalar tensor 
theory is the Brans-Dicke one which aimed to satisfy Mach ideas. 
Since the eighties some new justifications appeared mainly related 
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to particle physics theories. As instance, supersymmetry supposes 
the equahty between fermionic and bosonic degrees of freedom and 
needs several scalar fields to exist. The Higgs field is a scalar field. 
Some other reasons to believe that scalar fields exist are related to the 
cosmology: they could explain the flattening of some spiral galaxies 
rotation curves |2|( dark matter), the late time Universe accelerat- 
ing expansion^ E](dark energy) or the inflation. 
In this work we will be interested in a massive scalar field, non mini- 
mally coupled to the curvature. In p], it is called the Hyperextended 
Scalar Tensor (HST). It is such that the gravitationnal constant varies 
with respect to the scalar field. We will perform the calculations in 
the Einstein frame for which, after a conformal transformation of the 
metric, the field is cast into a minimally coupled and massive one. 
However it is then coupled to the perfect fluid. Consequently the 
matter does not follow the spacetime geodesies. 
The geometrical framework will be the homogeneous and spatially 
flat Bianchi type / model. It is an anisotropic model which gener- 
alises the flat Priedmann-Lemaitre-Robertson- Walker (FLRW) one. 
The Bianchi models allows understanding how the isotropisation ap- 
peared. If the isotropy and homogeneity of our Universe is well es- 
tablished until the decoupling periodfSj, one has to remember that it 
is a hypothesis concerning the early Universe. Another justification 
for an initially anisotropic state of the Universe is that the FLRW 
singularity approach is not generic. A generic approach could be 
oscillating as the one of the Bianchi type IX model. It has been 
conjectured by Belinskij, Khalatnikov and Lifchitz (BKL)[3IH] that 
it should be shared by the most general anisotropic and inhomoge- 
neous models. This conjecture has been revisited recently by Uggla 
and others)^. Here, we will study the Bianchi type / model whose 
singularity is not oscillatory but which is a spatially flat model in 
agreement with WMAP data. 

Our goal will be to study the isotropisation process of the Bianchi 
type / model in presence of a non minimally coupled and mas- 
sive scalar field. To this end, we will use the ADM Hamiltonian 
formalism ^01 It allows writing the field equations as a first order 
differential system. We will use the dynamical systems methods|12j 
to find its isotropic equilibrium states. The plan of the paper is as 
follows. In section 121 we write the Hamiltonian field equations. In 
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section IHl some necessary conditions for isotropisation and the be- 
haviours of the metric and potential when such a state is reached are 
described. In the last section, we summarize our results and study 
the isotropisation of the Brans-Dicke and low energy string theories 
when the potential is a power or an exponential law of the scalar 
field. 



2 Field equations 

The action of the HST in the Einstein frame writes: 

S = (16^)-i J[R- (3/2 + oj)cl)'^'cl)^^(l)'^ - U] V^d^x + </>) 

(1) 

4> is the scalar field, u and U are respectively the Brans-Dicke cou- 
pling function and the potential. They depend on (p. Sm is the action 
standing for a perfect fluid coupled to the scalar field. Its equation 
of state is Pm = (7 — l)Pm with 7 E [1,2]. The metric for the Bianchi 
type I model is: 

ds^ = -df + RlgijUj'J (2) 

The gij are the metric functions and tOi the 1-forms specifying the 
Bianchi type I model. In order to use the Hamiltonian formalism, 
we rewrite this metric with a 3+1 decomposition of spacetime: 

ds^ = -{N^ - NiN')dn^ + 2NidnJ + Rle-'^^^'^^'^JJ (3) 

N and A', are the lapse and shift functions. describes the isotropic 
part of the metric and will be considered as a time coordinate. We 
will show latter that it is a monotonic function of the proper time t. 
Moreover, we define a comoving 3- volume V as V = e~^^ . The f3ij 
stand for the anisotropic part of the metric, they were parameterised 
by Misner^^ in the following way: 

= diag{(i+ + \/3/?_ , 13+ - V2,f3- , -2/?+) (4) 

pi = 27r4 - 2/37:614 (5) 
6pij = diagip+ + V3p-,p+ - V^p-, -2p+) (6) 
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The pij are the f5ij conjugate momenta. To find the ADM Hamilto- 
nian, we rewrite the action as: 

S = (levr)-! / {tt'^^ + tt'I^^ + n^^-NC'- N,C')d*x (7) 
^ ' J ^ dt dt dt ' ^ ^ 

7r0 is the scalar field conjugate momentum. N and Ni are similar 
to some Lagrange multipliers. By varying the action with respect to 
these quantities, we find the constraints Co = and = with: 



(3)^ 2 

6 and A are respectively a positive constant and a scalar field function 
describing the coupling between the scalar field and the perfect fluid. 
The action Q may be derived from the one of a non minimally cou- 
pled scalar field as described in the appendice. Then the gravitation 
function G{(j)) is variable and we have the relation A oc ^^(^-s-t) . The 
energy conservation of the perfect fluid writes pm = XV . Hence, 
we will assume that A is a positive function of cp. The constraint 
Cj = is identically satisfied whereas the constraint Co = gives 
the ADM Hamiltonian: 

H^=pI+ pi + 12^^ + 247r2i?6e-6^C/ + ^6^(^-2)^^ (8) 
3 + loj 

When A = const, we recover the Hamiltonian when the scalar field 
is not coupled to the perfect fluid. For more details on Hamiltonian 
formalism, see [12]. The Hamiltonian equations then write: 

. dH I2(t?p^ 
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(12) 

A dot means a derivative with respect to O and a subscript cp a 
derivative with respect to the scalar field. We rewrite these equations 
with some bounded variables. For this we define: 

X = H-^ (14) 

y = e-^^VUH-^ (15) 

z = p<^0(3 + 2w)-i/2^-i (16) 

These new variables will be real if C/ > and 3 + 2ci; > 0. Thus the 
scalar field respects the weak energy condition. Each variable have a 
physical interpretation: 

• is proportional to the shear parameter S defined in 

• is proportional to {p^ — P(j,)/{dQ/dt)'^ , {dQ/dt)'^ being the 
Hubble constant when the Universe is isotropic, p,^ and the 
density and pressure of the scalar field. 

• z'^ is proportional to (p^ + p^) / {dil/ dt)"^ . 

• Prom the two last points it comes that the density parameter 
of the scalar field, Vt^ oc / (d^l / dt)'^ , is a linear combination 
of and z^. When the scalar field is quintessent these two 
variables are proportional to i}^. 

From the equation ((HI), we get: 

p^x^ + i?2y2 ^ ^2^2 + A;2 = 1 (17) 

where we put to simplify 

and we define the constants p'^ = p\ + p'^_ and = I^t^'^Rq. The 
equation (|T7|l may be considered as a constraint equation. It shows 
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that the variables x, y, z and k are bounded, k is not a new inde- 
pendent variable but is related to x, y and this last variable being 
able to diverge. It is proportional to the perfect fluid density param- 
oc pm/ {dO-/ dt)'^ and it may be rewritten under the following 
useful forms: 

fc2 = 5y'^U~^XV-'^ (20) 
Using the variables H14II16|) . the field equations become: 

x = m^y'^x-2,/2{-i -2)k'^x (21) 

y = y{Qiz + 2,B?y^ - 3) - 3/2(7 - 2)k^y (22) 
p 

z = i?2y2(32 ~ 2^ ~ ^'^^'^^ ~ ^''^^^ ~ l/2C/=^ (23) 

where the quantities £ and £m are defined by ^ = 4>U^U~^ {3+2ij)~^/'^ 
and £m = i;^A(^A~^(3 + 2a;)~^/^. £ and Im look each other because 
of the similar roles of U and A in the Hamiltonian Q. Both are 
multiplied by an exponential of The equation for (j) will be written 
as: 

(/) = 12z- ^— - (24) 

(3 + 2a;)V2 ^ ^ 

Summarising, the seven equations of the Hamiltonian system ipiT^ 
are reduced to a system of four equations (|21I24|I . It describes the 
evolution of four variables of which three are bounded. It comes ow- 
ing to the fact that, for the Bianchi type / model, the hamiltonian 
equations immediately give p± = const implying /5+ oc More- 
over, we will choose a diagonal form for the metric, i.e. Ni = 0. It 
allows getting N = 12iTRl,H-^e-^^ with dt = -Ndn. 

3 Stable isotropic states 
3.1 Defining isotropy 

Following Collins and Hawking|14|. isotropy arises when O — > — oo. 
Universe is thus forever expanding Universe with df3±/dt oc e'^^ 0. 
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Moreover, defining aij = {de^^ / dt)m{e~^) j-^^ aiid a"^ = aijaij, we 
must have ^^^^ — > 0. This last condition means that the anisotropy 
measured locally through the constant of Hubble tends to zero. It 
implies that the shear parameter x cc P± = dj3± / dtdt / dO. 0. Con- 
sequently, isotropisation occurs when x vanishes in O — > —oo. It is 
thus a stable state taking place for a diverging value of t. 

We deduce that the Universe can become isotropic following three 
different ways. We name them respectively class 1, 2 and 3. They 
are defined as follows; 

• Class 1: in the vicinity of the isotropy, all the variables (x, y, z) 
reach equilibrium with y ^ 0. It is generally possible to de- 
termine the asymptotical forms of the metric functions and 
potential whatever the forms of uj, U and A. 

• Class 2: in the vicinity of the isotropy, all the variables (x, y, z) 
reach equilibrium with y = and thus /c^ — > 1 — 122:^. Until 
now, we did not succeed in finding the Universe asymptotical 
state. A numerical example will be shown in the last section. 

• Class 3: in the vicinity of the isotropy, only the variable x 
reaches equilibrium but not necessarily the others variables. If 
y and z do not reach equilibrium when O — > — oo whereas they 
are bounded, they necessarily oscillate without damping. Thus 
their first derivatives oscillate around zero. This phenomenon 
will arise if I or /and tm sufficiently oscillate when — > — oo 
such that the signs of y or /and z change continuously. Of course 
in this case, once again it seems difficult to determine the Uni- 
verse asymptotical state (but we hope not impossible). Class 
3 isotropisation has been observed numerically in presence of a 
complex scalar field in jl^. 

In this paper, we will consider the first class which is agreement with 
the observations. Moreover, itallows a full description of the Universe 
asymptotic state when it isotropises. 
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3.2 Assumptions for the stability of our results 

In this paper, we will determine the isotropic equilibrium points, 
some necessary conditions for isotropisation and the asymptotical 
behaviours of some functions in the neighbourhood of these points. 
However we will get these behaviours by neglecting two quantities 
in the vicinity of the equilibrium. The first one is the variation of a 
function / of the scalar field whose form depends on the asymptotical 
behaviour of k. The second one is the variations of the variables 
(y, z, k). In other words, we will assume that all these quantities tend 
sufficiently fast to their equilibrium values. We had already talked 
about this problem in and we reproduce below the discussion of 
this last paper. 

The first type of assumption is related to £ and £m- Let f{£,irn) be 
a function of the scalar field, that we will define below, and tending 
to a constant /o in the vicinity of the isotropic state. We will assume 
that this function reaches sufficiently quickly its equilibrium value 
/o, i-e. 

• When / tends to its constant equilibrium value /o (vanishing 
or not) such as / — > /o + 5f, J (/o + 5f)dQ — > fo^ + const. 

We will check this assumption each time we will use our results. If 
it is not true, the asymptotical behaviours for the metric functions 
(and potential) will be different from the laws we will derive below. 
This problem could be overcame since our results allow to calculate 
4>{0,) and thus f{0,). Hence, it should be easy to generalise them 
by keeping the / fdO, term instead of considering that it tends to 
foil, + const. However, they would not be on a closed form. 

The second type of assumption can not be raised so easily. In the 
same way, the asymptotical behaviours we will determine will be true 
only if the variables {y, z, k) tend sufficiently fast to their equilibrium 
values. It means that we have to make the same kind of assumption 
for (y, 2, k) as for /. A solution to solve this problem would be to 
consider some small perturbations of these variables in the vicinity 
of the equilibrium. However, until now we did not succeed to get any 
interesting results, even for an empty flat model. 
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To summarize, the results of this paper concerning the asymptotical 
behaviours will be valid for a class 1 isotropisation if the functions / 
defined below and the variables (y, z, k) tend sufficiently fast to their 
equilibrium values. From a physical point of view, it means that the 
Universe tends sufficiently fast to its isotropic state. The assumption 
on / may be easily raised but the ones on (y, z, k) need a more careful 
examination. 



3.3 Asymptotic state when k 

In what follows, we look for the isotropic states with /c / 0. The case 
with A; — > will be analysed in the section [^TKl 

Equilibrium points 
First we calculate the equilibrium points of the equations system l(2Tl 
1^ . Then, we introduce them in the constraint lfT7|l to find k. The 
equilibrium points write: 

- 1) - 97(7 - 2)] / , md^) ' 

^2 = (0,±— -[4^„(^„-^)-37(7-2)]i/2 7 ^ 



Eq belongs to the class 2 and is thus discarded. For the two other 
points, we find for k: 

^ 2£{£ - Im) - 37 



/c is a real as long as: 

^(^ - ^m) > \l (25) 

This inequality implies a real value for points Ei which are thus 
discarded from further considerations. Consequently, the only equi- 
librium points corresponding to an isotropic class 1 stable state are 
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Figure 1: Equilibrium point E2 in {x,y,z) = (0,0.81,-0.08) when 
A: / and {R, 7, £, Im) = (1, 1, -1-23, 1.58). 



the E2 ones. A numerical simulation representing one of them is 
represented on figure Q They are real and bounded if respectively: 

4C(C-^) > 37(7-2) (26) 

i -h (27) 

i.e. [/ 7^ A. The first condition is always satisfied when there is no 
coupling between the matter and the scalar field(^m = 0). We will 
show in section 13.41 that t and tm can not diverge but at the same 
order and that k tends to a non vanishing constant in the vicinity of 
the E2 points. 

Monotonic functions 
The equation ipijl shows that x is a monotonic function of constant 
sign. We deduce that the metric functions, whose derivatives with 
respect to the proper time express as some linear functions of x, can 
only have a single extremum. Prom the form of the lapse function N 
and the relation dt = —NdQ, it comes that O is a monotonic func- 
tion of the proper time. Its value in —00 corresponds to late time 
epoch when the Hamiltonian is initially positive. We will find the 
same monotonic functions when A; — > 0. 
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Asymptotical behaviours 
The function / (see subsection l,S.2p is defined in this subsection as 
/ = £(i — £m)~^ ■ Then hnearising the equation pT]! for x in the 
vicinity of the isotropic equihbrium state, we find: 

3[2^m+^(7-2)l o 

X XqC 2{f-£„) ^ 

with xq an integration constant^, x vanishes when — > — oo if 
the reahty conditions for k and E2 are respected. Prom the lapse 
function N and the relation dt = —NdQ, we find that tends to 
the increasing function of the proper time: 

e-^i ^ t (28) 

when 2(e-'e ) tends to a non vanishing constant. This is always true 
since, as it will be shown in the next section, £ and im can not diverge 
but at the same order with i ^ Im- Moreover, the reality condition 
l(25]l for k shows that I does not vanish. Since in the next section 
we will also prove that y can not vanish, from y definition we derive 
that the potential asymptotically disappears as ?7 — > 
To check the necessary conditions for isotropy H25B27|1 when a;(i;A) and 
[/((/>) are specified, we need to know the asymptotical behaviour of 
0. For this, we write (|2^ in the neighbourhood of the equilibrium. 
Then, we get that asymptotically 4> behaves as the asymptotical so- 
lution of the differential equation: 

'^ = 37(^-^)-^ (29) 

when Q — > —00. This asymptotical equation is a good illustration of 
the assumptions on (y, z, k) described in subsection l,S.2> to determine 
l(29|l . we have replaced z in (|2^ by its equilibrium value jjjz^-^^ 
neglecting any small variation bz when z approaches equilibrium. 
However, if it tends to ^^^J^ ^ slower than then bz should be 
taken into account in (pi^l instead of being neglected. 



3.4 Some important results for the points 

Integrating near equilibrium leads to 

^This result is always valid if I and tm diverge at the same order. 
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Uq being an integration constant. We then deduce from the expres- 
sion IpUl that k tends to a non vanishing constant as long as y / 
what is always true. Indeed, the only ways for y to vanish are if 
t » Irri ox Mm{(-m — ^) ^ '^lil — 2). In the first case, the con- 
straint or equivalently l(28)l shows that /c — > 1. But if we consider 
the form for k and the limit U UqXV^^ , it comes that if 
y^0=^A;^07^1. For the same reason, y can not tend to when 
4^m(^m — ^) ^ 37(7 — 2). It follows that y is never vanishing in the 
vicinity of E2 . We get a similar result if we consider the divergence of 
tm when tni » L In this case, y tends to a non vanishing constant 
and A; to 0, which is in disagreement with the form (PT)]! of k and 
the limit U — > Uq\V~'^ . On the other hand, if £ and Im diverge at 
the same order without converging one to the other, y and k tend to 
some non vanishing constants and the constraint is respected. Since 
A oc UV^ and y is not vanishing, we deduce from lfTK|l that 

A ^ e t-^^ 

and from (jSHl that 

A t~^~ 

Thus from the reality condition for A; , it comes A > t ^ aj^-' 
In the same time, the energy densities for the perfect fluid and 
scalar field write respectively p = \V~"' U and = |7^(^~^ — 
+ ^U. When £ and im tend to some constants, since 
U P(j, (X p^ (X pm- the scalar field and the perfect fluid 

energy densities behave in the same way. If both £ and £m diverge at 
the same order, the kinetic term in p^p is larger than p,p » Pm 
and the energy density of the scalar field dominates. 

3.5 Equilibrium point when A; ^ 

This section is divided in three parts depending on £m = strictly, 
£mk'^ ^ or £mk'^ 7^ 0. 



In this first part, we recall and complete the results got in ^21 when 
no coupling exists between the scalar field and the perfect fluid. 
When there is no perfet fluid, fc = strictly and the reality condition 
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for the equilibrium points writes £^ < 3. Near isotropy, the metric 
functions tend to when £ tends to a non vanishing constant or to 
an exponential when I vanishes. When /c 74 0, the reality condition 
for the equilibrium points is > 3/27 and the metric functions tend 

to t^T . When A; — > 0, we recover the same values for the equilibrium 
points as when no perfect fluid is present[THj but now, A; — > implies 
l"^ < 3/27. The asymptotical behaviour of the metric functions is 
the same as without a perfect fluid. We had not noticed this last 
inequality in ^2) nor that U — > when k — > const / 0. 
Once again, these results rest on the assumptions of subsection 13.21 
with now f = i'^ when /c = or /c — > 0. When they are not true, 
meaning that the Universe does not reach its isotropic equilibrium 
state sufficiently quickly, the asymptotical behaviours of U and 
are generally different. As instance when k = strictly and if £^ 
vanishes as nQ~^ with n < 0, the integral of / does not tend to 
a constant. Then we can show that the potential will diverge as 
(— and the metric functions will tend to exp j^^^t^^^'^'^^^ 
with n £ ] — 1,0[ such as the Universe be expanding. This solution 
is different from the classical solutions found when we neglect the 
variation of i near equilibrium. It shows that the assumptions on / 
that we will also use in the next sections, have to be checked each 
time we apply our results to a specific scalar-tensor theory. 

If fc — > with Imk'^ 0, again we recover the same equilibrium 
points and behaviour for x as in the absence of a perfect fluid. That 
is X ^ xoe^^~^^^^ with < 3 such that x ^ in — > —00 and 
the equilibrium points are real. Consequently, using the form ^ 
for k'^, it comes fc^ Xf^'^i^fij-^^)^ when Q —00. When £^ tends 
to a non vanishing (vanishing) constant smaller than 3, — > 
(respectively — > g^^^'^'^o^o)"^*). Hence, A; — > if 

Ae2(3/27-^^)n ^ (30) 



and thus A < t 2i^^ (respectively A < e3'i'(i2'^«o^o) ^ty ^^le 
same way, £mk'^ if: 
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Contrary to the case im = 0, the condition k —>■ does not automat- 
ically restrict the set of i allowing isotropisation: it is the form of 
A which will lay down the law. Since we consider a class 1 isotropi- 
sation such as y 7^ and k ^ 0, we have XV~'^ « U and thus 
P4> ~ P(j) » Pm- the scalar field energy density dominates the Uni- 
verse. 

The asymptotical behaviour of the scalar field when — > — cxd is 
given hy|lS|: 

■ ^ (31) 



C/(3 + 2uj) 

^k^7^0 



Since A; — > whereas (mk -f^ 0, it means that Im have to diverge. 
The equilibrium points when A write: 

^3 = (0,±i?-\0) 

They are such as k^ = —£i^. The dynamics approach these points 
in the same way as on the figure ^ k will vanish and will be real 
if respectively i « £m and < 0. Moreover, imk'^ 7^ if 

£ tends to a non vanishing constant or diverges with zi bounded. 
Mathematically, the £"3 point could be the asymptotical limit of E2 
when £rn diverges and I « £m- However, this divergence is forbidden 
by the constraint. 

Near £"3, we find that x e^^ , showing that the Universe tends to 
a De Sitter model, i.e. g(i27ri?gxo)t^ ^.j^g potential to the 

constant {Rxq)~'^ . As previously, we have then ^ if 

Ae^T^ ^ 

i.e. A < e^'^'^^^'^^o^o) In the same way, Imk"^ does not vanish if: 

£„Ae3^^ /> 

Again, y being different from and considering the form (|2()|) for 
A;^, we have pm = XV « U such as A; — > and thus pm « 
P<j) ~ V4>'- the Universe is scalar field dominated. From (|T9|l and the 
limit of k near equilibrium, we determine the scalar field asymptotical 
behaviour: 

5J_^ = e3^" (32) 
A^ U 
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4 Discussion 

The discussion is divided in three parts. In the first one we summarize 
our results and in the second one we consider some apphcations for 
Brans-Dicke and low energy string theories. We conclude in the third 
one. 

4.1 Summary 

We have studied the necessary conditions which may lead the Uni- 
verse to a class 1 isotropisation. It depends if k does not vanisli, 
vanishes with £mk'^ — ^ or with Imk"^ -h 0- We assumed that 3 + 2a; 
and \J are some positive functions of the scalar field and that the 
isotropic state was reached sufficiently fastly. Below we summarize 
our results. 

Case. 1: k-Z^O 

We define the quantities £ = (l)U(f)U~^{3+2uj)^^/^ andi,n = (i>X0X^^(3+ 
2uS)~^l'^ . Letptj), Pff) and pm he respectively the pressure and density of 
the scalar field, the density of the perfect fluid. Some necessary condi- 
tions for Bianchi type I isotropisation in presence of a massive scalar 
field minimally coupled to the curvature hut not minimally coupled to 
the perfect fluid are: 

• i trn (equilihrium points are hounded) 

• 4:im{£m — > 3(7 — 2)7 (reality condition) 

• £{£ — £rn) > §7 (reality condition) 

• £ and Im are hounded or diverge in the same way (the constraint 
is respected) 

When isotropy is approached, the metric functions behave as t ^ 
A — > f'^'f' whereas the potential decreases as t^^ . When £ and £m 
do not diverge, there is an equilihrium hetween the scalar field and 
the perfect fluid: p^f, oa p^j) oc pm- When both diverge, Ptp — P(j> » Pm 
and the Universe is scalar fi,eld dominated. Asymptotically, the scalar 
field checks the relation U —>■ UoXe^'^^ . 
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This last relation determines the asymptotical form of (j) and thus 
these of I and im- Note that in the case = 0^2) the metric 
functions asymptotical behaviour does not depend on 4> and is al- 

2 

ways t^i . It thus prevents any late time acceleration. Hence, it is 
the existence of a coupling between the scalar field and the perfect 
fluid which allows the appearance of an accelerated expansion when 
k — > const / 0. Then, since when ^ and are bounded we have 
P<)) oc Pm, it follows that oc Qm and the coincidence problem could 
be solved. 

Case 2: k ^ and £mk'^ 

We define the quantities i = (j)U^U~^{3+2uj)~^^'^ andim = 4>Xtf,X~^{3+ 
1io)~^l'^ . Some necessary conditions for Bianchi type I isotropisation 
in presence of a massive scalar field minimally coupled to the curva- 
ture but not minimally coupled to the perfect fluid are: 

• < 3 (reality condition) 

• Ae2(3/27-^')f^ ^ (Condition for k ^ 0) 

• CAe(37-2^')f^ ^ (Condition for l^k'^ Q) 

If £^ tends to a non vanishing constant, the metric functions tend to 
t^ and the potential vanishes as If l"^ vanishes, the Universe 
tends to a De Sitter model and the potential to a constant. In any 
cases P(j)—P(j) >> Pm o,nd the Universe is scalar field dominated. The 
asymptotical behaviour for the scalar field is this of the asymptotical 

solution of<j) = 2 ^^3^^) . 

These results include the ones got when there is no perfect fluid |18|. 
For sake of clarity, we chose to express the above (as well as below) 
limits for A; ^ and £rnk ^ as some functions of and (/>. How- 
ever, these two last quantities are asymptotically defined with respect 
to the proper time t by the behaviours of the metric functions and 
potential. 

Case 3: k ^ and imk'^ -/^ 

We define the quantities i = (j)U^U~^{3+2uj)~^^'^ andim = 4'X^X~^{3+ 
2u;)^^/^. Some necessary conditions for Bianchi type I isotropisation 
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in presence of a massive scalar field minimally coupled to the curva- 
ture but not minimally coupled to the perfect fluid are: 

• Im diverges and t — > const ^ or diverges such that z£ ^ 
(condition for imk"^ ^ ^) 

• £ << £m or Xe^^^ (condition for k ^ 0) 

• ££^ < (reality condition) 

The Universe tends to a De Sitter model and the potential to a con- 
stant. Since P(f>—P(f> » Pm, the scalar field asymptotically dominates 
the Universe and checks the equation ^j^^ = e'^'*^^- 

The cases with k -/^ Q and A; — > are strictly separated by the asymp- 
totical behaviour of A since the first one impHes A > t ^ ^j^-^ g^j^^j 
the second one A < r^^^"^^^^ (or A < eS-yCis^rflg^o)-!* ^j^g^ £ 0). 
The two cases with /c ^ are distinguished by the fact that for the 
first one, the Universe tends to a De Sitter model when £ — > and 
the second one when £ ^Q. 

4.2 Applications 

In what follows, we apply our results to some Brans-Dicke and low 
energy string actions. For that, we use a conformal transformation of 
the metric described in the appendice. It casts the minimally coupled 
scalar-tensor theory Q in the Einstein frame where our results take 
place, into a non minimally coupled scalar-tensor theory H,S,S|1 in the 
Brans-Dicke frame. Obviously, when isotropy arises in the Einstein 
frame, it also occurs in the Brans-Dicke frame. Thus the necessary 
conditions for isotropy are the same in both frames. However, the 
metric functions generally behave differently. 

We will illustrate each of the four applications below with some fig- 
ures representing the behaviours of x, y, z, /c, (f) and £ in the Einstein 
frame and in the Q. time. The numerical integrations will be done 
using the initial conditions = 0.14, yo = 0.25, zq = 0.12. xq is 
calculated with the constraint (|T7|l where we put p\ + p'^_ = = 1, 
R = 1 and 6 = 1 (the constant in the definition of A;). The behaviours 
of the metric functions a, (3 and 7 in the Brans-Dicke frame, and their 
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derivatives, will be also represented but in the proper time t with ini- 
tial conditions oq = —1.53, /5o = —1-25, 70 = 0.12, dao/dro = 2.48, 
dPo/dTQ = 1.55 and dj/dro = 0.33. Here, the r time is defined 
as dt = Vdr. To get the figures, we numerically integrated the 
field equations of the Lagrangian formalism using a 5 order Runge- 
Kutta method^. d4>Q/dTQ is calculated using the constraint equation 
of this formalism. Each time, a dust fluid and a null initial time were 
asumed. 



Brans-Dicke theory with an exponential potential 

Consider the class of theories defined by (0) such that: 



U 
A 



Using the conformal transformation, it can be cast into the non min- 
imally coupled scalar field theory (|33|l defined by: 



G 



U 



(;/)3{4-37) 

5(1 



m 



-2(H 



9(4 - 37)2 

) „n</> 



03(4-37) 



3(4-37) 



The Brans-Dicke theory with an exponential potential is recovered 
for m = 3(37 — 4). 

The quantities £ and Im are defined by: 



ncf) 



V3 + 2^0 
m 



V3 + 2^0 

£m can not diverge and consequently the case 3 never happens. More- 
over 3 -|- 2ojQ have to be positive. For the case 1 . near the isotropic 
equilibrium state, we have for the scalar field: 



^n4>i-(2+m) 



37Q 



^We used Java oriented object programing to perform the calculations. 
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Since i is bounded, (/> can not diverge and should asymptotically 
vanish. It implies that m < —2 and finally (j) fpjjg 
second reality condition gives: 

4(2 + m) - 37(3 + 2loo) 

2(3 + 2wo) ^ 
But since m < —2, 7 > and 3 + 2ujo > 0, it can not be satisfied. 
Consequently, a class 1 isotropisation does not arise. 

Let us consider the case 2 . Integrating the differential equation for 
(f), we get: 

4> -in- 

n — (/)oe^+^"o 

Then, when Q — > —00, (p 2n~^, £ — > and A tends to the constant 
(2^^^)"*, implying n > 0. If the Universe isotropises, it will tend to 
a De Sitter one. Remark that (p, and thus n, have to be positive such 
that A be a real function. 

Using the conformal transformation, when isotropisation occurs in 
the Brans-Dicke frame where (f) is non minimally coupled to the cur- 
vature and since A tends to a constant, the metric also tends to a De 
Sitter form (see figure [SJ- 

A class 2 isotropisation is also possible when n < and is plot- 
ted on figure 01 As above noted, such a range for n is impossible 
for a class 1 isotropisation since A would be a complex function. It 
is the only example of class 2 isotropisation we have found until now. 

Brans-Dicke theory with a power potential 

Consider the class of theories defined by the Lagrangian Q such 
that: 

u = 4)'^ 
\ = (T 

If we apply again the conformal transformation, we obtain the non 
minimally coupled scalar tensor theory defined by: 

m 

G = (/)3(4-37) 
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Figure 2: These figures represent the approach for a class 1 isotropisation when 
uio = 2.3, n = 1.5 and m = 1.1. As expected, x tends to 0, (p to the constant 2/n = 1.33 
and £ (here named ell) to 0. The convergence of </> to a constant is in accordance with 
the fact that U also tends to a constant and the Universe to a De Sitter model. In the 
Brans-Dicke frame, the derivatives of the metric functions a, f3 and 7 tend to the same 
behaviour: isotropisation occurs. 




Figure 3: These figures represent the approach for a class 2 isotropisation when 
ojo = 2.3, n = —3.1 and rn = 1.1. x always tends to but also y. 4> and thus i diverge. 
Note that ij>, y, z and £ undergo damped oscillations. In the Brans-Dicke frame, the 
derivatives of the metric functions tend to the same behaviour showing isotropisation. 



LO = 



3 _ 

r ~ 9(4-37)2- 



Wo J 



03(4-37)' 



u 



(f) 3(4-37) 



The Brans-Dicke theory with a power law potential is recovered for 
m = 3(37 - 4). 
We calculate that: 

e 



n 



V3 + 2loo 

m 

V3 + 2^0 
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with 3 + 2uJo > 0. Anew im can not diverge and the case 3 is ex- 
cluded. For the case 1 . it is necessary that n / m such that £ -/^ im- 
Asymptotically the scalar field behaves as: 

Consequently, in O ^ — oo, (p 0(0 diverges) if m — n < (respec- 
tively m — n > 0). The reality conditions write: 

4m{m - n) + 37(2 - 7)(3 + 2ujo) > 

2n(n - m) - 37(3 + 2ujo) > 

The second one will be respected if n>0(n<0) when cp (respectively 
(j) diverges). We find then that if an isotropic state is reached, the 

2(n-m) 2m 

metric functions tend to t and A to t " . 

Using the conformal transformation, we deduce for the non mini- 
mally coupled theory that the metric functions will tend to: 

m(8-57)+2n(37-4) 
I 7[m,+3n(37-4)J 



All these behaviours are illustrated on figure 01 




Figure 4: These figures represent the approach for a class 1 isotropisation when 
LOO = 2.3, n = —3.1 and m = 1.1. Anew x tends to vanish, y and k to some non 
vanishing constants showing that U oc Ae^''^. if> diverges since m — n > 0. In the 
Brans-Dicke frame, the Universe isotropises. 



For the case 2 . we get for (p: 



(p — > g 3+21^0 
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Hence k will vanish when Q — > — oo if 2n(m — n) + 3j{3 + 2iOo) > 0. 
The reality condition for the equilibrium points will be respected if 
n^(3 + 2wo)"^ < 3. The metric functions then tend to t{3+2a;o)n-2 
when n 7^ or to a De Sitter model when n = 0. 

In the Brans-Dicke frame where the scalar field is non minimally 
coupled to the curvature, the metric function will tend to: 

mn + 3(37-4)(3 + 2a;o) 
^ n|m+3n(37-4)] 



when n 7^ 0. If n = 0, the behaviour of the metric functions is the 
same as in the Einstein frame and the Universe tends to a De Sitter 
model. This case is illustrated on figure 




Figure 5: These figures represent the approach for a class 1 isotropisation when 
LOO = 2.3, n = 1.5 and m = 1.1. Here, k tends to vanish and the scalar field energy 
density dominates the Universe. 



Low energy string theory with an exponential potential 

We consider the theory defined by (Q) and such that: 

U = e""^ 
X = e""'' 

Applying the conformal transformation, we define the following non 
minimally coupled scalar tensor theory: 

m I 
G = e 3(4-37) 
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2 18(4 - 37)2 



3(4-37)' 



The low energy string theory with an exponential potential is then 
recovered when m = 3(4 — 37), ujq = 5/2 and oji = —3/2. 
We calculate I and im and we obtain: 



V3 + 2(^.2^0 + 2a;i 
mcf) 



V3 + 2</.2c^o + 2wi 
These expressions show that we will never have £ « im and thus 
the case 3 never occurs. For the case 1 . it is necessary that m ^ n. 
Moreover, we find for the scalar field: 

37^1 



4>^ 4>o + 



n — m 

Hence, (j) diverges and £ and im tend to some constants. They will 
be real if ujq > 0. The reality conditions write: 

2m(m-n) + 3(2 - 7) > 

n{n — m) — 2>"fuJo > 

Wo being positive, the second condition needs n{n — m) > and thus 
n 7^ 0. Consequently, when isotropisation arises, the metric functions 

2 n — m _2 — 

and A respectively tend to t and t n . 

We deduce that in the Brans-Dicke frame, when isotropisation arises, 
the metric functions will tend to: 

m(8-57) + 2n(37-4) 
I 7[m + 3n(37-4)| 

This case is represented on the figure El 

Concerning the case 2 . the scalar field asymptotically behaves as: 
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Figure 6: These figures represent the approach for a class 1 isotropisation when 
UJO = 2.3, n = —3.1 and m = 1.1. k tends to a constant showing the equilibrium between 
the scalar field and the perfect fluid. Remark the existence before this equilibrium of 
a period during which the density of the scalar field dominated the one of the perfect 
fluid. 

Consequently, depending on the sign of the square root, we have two 
branches such that (/> — > or (/) — > nuj^^fl. For the first one, £ — > 
and the Universe tends to a De Sitter model. The limit allowing 
/c — > is always respected. For the second one, £ n(2a;o)~^/^ and 
thus, isotropisation needs > and n^(2a;o)^^ < 3. If n / 0, the 

metric functions tend to t~nT and the limit allowing — > is satisfied 
if < If n = 0, the Universe tends to a De Sitter model and the 
limit A; ^ is always satisfied. 

Again, in the Brans-Dicke frame, we deduce that when isotropisa- 
tion arises and the scalar field vanishes or n = 0, the metric functions 
tend to the same form as in the Einstein frame because A tends to a 
constant. When the scalar field diverges and n / 0, they tend to: 

n^(97-13)+3(77-8)t^0 
f n2(97-13) + 37t^0 

Low energy string theory with a power potential 

We now consider the minimally coupled Lagrangian Q with: 

LO = UJo<p^ + UJl 

U = (Ife""^ 
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Figure 7: These figures represent the approach for a class 1 isotropisation when 
Wo = 2.3, uJi = 0.5, n = 1.5 and m = 1.1. k tends to vanish showing that the scalar field 
energy density dominates the one of the perfect fluid. In the same way, cp vanishes. 



Applying the conformal transformation, it is cast into the following 
non minimally coupled theory: 



G = e 3(4-37) 

■3 

2 



LO = 



3m 



U 



.2 

)<f> 



18(4 - 37)2 



/,g 3(4-37) ^ 



The law energy string theory with a power potential is recovered 
when m = 3(4 — 37), n = 2, loq = 5/2 and loi = —3/2. 
Calculating £ and im, we get: 

p + ncf) 



V3 + 2(t>'^uo + 2wi 

^ ^ rn(^ 

v/3 + 2(/.2c^o + 2uJi 
Again, it is impossible that diverges and in the same time i « 
£m- Thus the case 3 is excluded. For the case 1 . we show that the 
scalar field behaves as: 

</, = p(m - n)-^ProductLog{{n - m)e^'^P~''^^-'^°^) 

When pj~^ > 0, the scalar field vanishes, otherwise it diverges. 
Then, (n — m)p~^ have to be positive otherwise is complex. 
When ^ — > 0, it is necessary that 3 + 2a;i > such that I and Em be 
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real and the reality conditions for the equilibrium points reduce to 

2 

2p^ — 37(3 + 2uJi) > 0. Then, the metric functions tend to and A 
to a constant. This case is plotted on figure |H1 
When (/> — > 00, it is necessary that cjq > such as £ and im be real 




Figure 8: These figures represent the approach for a class 1 isotropisation when 
LJo = 2.3, Lui = 0.5, n = —3.1, m = 1.1 and p = 3. k oscillates to a constant and (f> 
vanishes. Note the strong oscillations of y, z and k. 



and n ^ m such as i does not tend to 
the equilibrium points write then 



The reality conditions for 



2m{m — n) + 37(2 — 7)11^0 > 

n(n — m) — 37u;o > 
It implies that n(n — m) > and n ^ 0. The metric functions tend 

to t 3"T and X —>■ t n , Some figures similar to the figure |H1 but 
with diverging (j) may be obtained. 

In the Brans-Dicke frame, the metric functions tend to the same 
form as in the Einstein frame during isotropisation if — > 0. When 
(j) diverges, they tend to: 

m(8-57) + 2n(37-4) 
I 7[m + 3n(37-4)| 



Let US examine the case 2 . The scalar field is such that: 
"(3 + 2a;i)ln</) n'^{3 + 2uJi) + 2p'^iJo 



)+l/2 



P 



ln{p + ncf)) + 



2u;q(/) 



n 
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Hence, it exists three cases such that O — > — oo. 
In the first one, (j) tends to vanish and it is then necessary that 
59 > and 3 + 2cji > 0. £ ^ p(3 + 2u:i)-^/'^ and thus we need 
j92(3 + 2u;i)~i < 3. The metric functions tend to t{3+2'^i)/p\ k al- 
ways tends to as long as £^ < 3/27. This case is shown on figure El 
Since 4) vanishes, A tends to a constant and the results are the same 
in the Brans-Dicke frame. 
In the second one, (j) diverges as nS"^- niust be positive and 




Figure 9: These figures represent the approach for a class 1 isotropisation when 
u)Q = 2.3, uji = 0.5, n = —3.1, m = 1.1 and p = 0.7. k and vanish, t tends to 0.35 
which is smaller than 3/27 = 3/2 

Wo > thus implying (j) — > +00 and n < 0. Then, I tends to 
n{2ujQ)~^/'^ and it follows that a necessary condition for isotropisa- 

tion is v? {2ijjQ)~^ > 3. Then, the metric functions tend to and 
A: to if n{m — 2n) + Gjujq > 0. In the Brans-Dicke frame, we find 

that the metric functions tend to t 

In the third one, cf) — > —pn~^ which implies [— n^(3 + 2uJi) — 2k'^ujQ] (pn^)^^ > 
0. Then, £ — > and the Universe tends to a De Sitter model. The 
condition A; — > is always respected. Once again A tends to a con- 
stant and in the Brans-Dicke frame, the metric functions tend to the 
same form as in the Einstein frame. 

4.3 Conclusion 

We have found some necessary conditions for isotropisation of Bianchi 
type / model with a massive scalar field, minimally coupled to the 
curvature but not to the perfect fluid. They depend on the asymp- 
totical behaviours of k and the product kim- We have then deduced 
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the asymptotical behaviours of the metric functions and the potential 
in the vicinity of the isotropy. A possible solution to the coincidence 
problem has also been found. Through some applications, we have 
shown how to extend our results to a scalar-tensor theory. The nec- 
essary conditions for isotropisation are the same in the Einstein or 
Brans-Dicke frames but the asymptotical behaviour of the metric 
functions are different. They may be determined via a conformal 
transformation of the metric. We have thus studied the isotropisa- 
tion of the Brans-Dicke and low energy string theories with some 
power or exponential laws of the scalar field for the potential. 

Parts of the calculus and phase portrait diagrams have been made 
with help of the marvellous DynPack 10.69 package for Mathematica 



4 written by Alfred Clark I http:/ / www.me.rochester.edu/courses/ 



ME406/webdown/down.html for download). 

A Perfect fluid conservation law when it is 
non minimafly coupled to the scalar fleld 

In this appendice, we calculate the energy momentum conservation 
law of the perfect fluid when it is non minimally coupled to the scalar 
field. This calculus is also made in pHI and more particularly in |2n| . 
Let us consider the Lagrangian of a non minimally coupled scalar 
field also known as hyperextended scalar tensor theory [5]: 

L = {G-'R - uj4>-'(I),^cI)'^' -U + T^fSg^^)^ (33) 

Then we define a conformal transformation of the metric: 

9ai3 = Gg^p (34) 

dt = VCdt (35) 

The frame related to gap is usually called the Brans-Dicke frame 
whereas the one related to gap is called the Einstein frame. In both 
cases, t and t are the proper times such as the 00 metric functions 
components are —1. Applying the transformation (|35|l casts the La- 
grangian 

(ESI into: 

L=[r- {3/2{G-')lG^ + ojG(t>-')(t>,f,(f>^'' - G^U + G^T^^Sgap] 

(36) 
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where the scalar field is now coupled non minimally with the perfect 
fluid but minimally with the curvature. Consequently, it comes: 

f = G^T 

We deduce the following energy conservation law: 

= 3G,„G2r°^ (since T^J^ = 0) 

= SG^aG'^g'^^f 

dG 

T^J^ = -3^G~^T (since G = G{t)) 



Let us remark that in this law is interpreted as the action of 
a force on matter due to the variability of the rest masses. Conse- 
quently, matter does not follow the spacetime geodesies. To simplify 
the calculations, we put p* = G^p and p* = G^p. Hence, we have 
f"^^ = {p*+p*)u°'u^+g"'^p. Moreover, we have assumed^ = (7— 
Thus, it comes: 

p*Y'y _ ^3(4-37) 

From this last result and the expression for the Lagrangian Lm for 
a perfect fluid calculated in pages 48-52], we are able to deter- 
mine the form of Hm, the term describing the matter in the ADM 
Hamiltonian. Indeed, we have: 



= -87rRlNe-^^p 



-SirRlNe-^^p* 
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and consequently: 

Hm = -2A7r^ g^/^L„, = 1927:^ RIG^^^-^"'^ e^'-"''^^^ > (37) 
We will write symbolically this relation under the form Hm = 5X{(l))e^^"' 

References 

[1] Tonatiuh Matos, F. Siddhartha Guzman, and Dario Nunez. 
Spherical scalar field halo in galaxies. Phys. Rev., D62:061301, 
2000. 

[2] S. Fay. Scalar fields properties for flat galactic rotation curves. 
Astronomy and Astrophysics, 413:799, 2004. 

[3] S. Perlmutter et al. Measurements of Q and A from 42 Hight- 
Redshift Supernovae. Astrophysical Journal, 517:565-586, 1999. 

[4] Adam Riess et al. Observational evidence from Supernovae for 
an accelerating Universe and a cosmological constant. Astro- 
physical Journal, 116:1009, 1998. 

[5] Diego F. Torres and Hector Vucetich. Hyperextended scalar- 
tensor gravity. Phys. Rev., D54:7373-7377, 1996. 

[6] D.N. Spergel et al. First year wilkinson microwave anisotropy 
probe (wmap) observations: Determination of cosmological pa- 
rameters. Astrophys. J., 148:175, 2003. 

[7] V. A. Behnskii, I. M. Khalatnikov, and E. M. Lifshitz. A gen- 
eral solution of the Einstein equations with a time singularity. 
Advances in Physics, 31, 6:639-667, 1982. 

[8] V. A. Behnskii, I. M. Khalatnikov, and E. M. Lifshitz. Oscil- 
latory approach to a singularity in the relativistic cosmology. 
Advances in Physics, 19:525-573, 1970. 

[9] Claes Uggla, Henk van Elst, John Wainwright, and George 
F. R. Ellis. The past attr actor in inhomogeneous cosmology. 
\gr-qc/0304.002 , submitted for publication to Physical Review D, 
2003. 



30 



[10] Hidekazu Nariai. Hamiltonian approach to the dynamics of Ex- 
panding Homogeneous Universe in the Brans-Dicke cosmology. 
Prog, of Theo. Phys., 47,6:1824, 1972. 

[11] C. W. Misner. Phys. Rev., 125:2163, 1962. 

[12] J. Wainwright and G.F.R. Ellis, editors. Dynamical Systems in 
Cosmology. Cambridge University Press, 1997. 

[13] Michael P. Ryan. Hamiltonian cosmology. Springer- Verlag, 1972. 

[14] C. B. Collins and S. W. Hawking. Why is the universe isotropic. 
Astrophys. J., 180:317-334, 1973. 

[15] S. Fay and J. P. Luminet. Isotropisation of flat homogeneous cos- 
mologies in presence of minimally coupled massive scalar fields 
with a perfect fiuid. Class. Quantum Grav., 21:1849-1878, 2004. 

[16] S. Fay. Isotropisation of Bianchi class A models with a minimally 
coupled scalar field and a perfect fiuid. Class. Quantum Grav., 
21, 6:1609-1621, 2004. 

[17] S. Fay. Isotropisation of the minimally coupled scalar-tensor 
theory with a massive scalar field and a perfect fluid in the 
Bianchi type I model. Class. Quantum Grav, 19, 2:269-278, 
2002. 

[18] S. Fay. Isotropisation of Generalised-Scalar Tensor theory plus a 
massive scalar field in the Bianchi type I model. Class. Quantum 
Grav, 18:2887-2894, 2001. 

[19] T. Damour and K. Nordtvedt. Tensor-scalar cosmological mod- 
els and their relaxation toward general relativity. Phys. Rev., 
D48:3436, 1993. 

[20] Jacob D. Bekenstein. Are particle rest masses variable? theory 
and constraints from solar system experiments. Phys. Rev., D15, 
6:1458, 1977. 



31 



